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CO ' Abstract. We survey recent progress in the problem of recovering a tensor field 

. from its integrals along geodesies. We also propose several open problems. 
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O ■ 1. Introduction 

cn 

This paper surveys recent results on the integral geometry problem of recovering a 
tensor field from its integrals along geodesies. The most basic example of the kinds of 
' transforms studied in this paper is the X-ray (or Radon) transform in the plane, which 
■ encodes the integrals of a function / in over straight lines: 

roo 

Rf{s,u)= f{su + tu^)dt, seR,ueS\ 

J — oo 

Here is the rotation of u by 90 degrees counterclockwise. The properties of this 
transform are classical and well studied [20j. The X-ray transform forms the basis for 
many imaging methods such as CT and PET in medical imaging. 

A number of imaging methods involve generalizations of this transform. In seismic 
and ultrasound imaging one encounters ray transforms where the measurements are 
given by integrals over more general families of curves, often modeled as the geodesies 
of a Riemannian metric. Moreover, integrals of vector fields or other tensor fields 
instead of just integrals of functions over geodesies may arise, and these transforms 
are also useful in rigidity questions in differential geometry. We will give more specific 
examples after having defined the relevant transforms precisely. 
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The geodesic ray transform acts on tensor fields on a compact, oriented Riemannian 
manifold (M, g) with boundary of dimension dim (M) = n > 2. We denote by ( ■ , ■ ) 
the ^f-inner product of tangent vectors or other tensors, and by | ■ | the ^f-norm. Let u 
denote the unit outer normal to dM. We denote by SM — )■ M the unit-sphere bundle 
over M: 

SM=\J S,, S^ = {ve T^M : \v\^ = 1}. 

The set SM is a (2n — l)-dimensional compact manifold with boundary, which can be 
written as the union d{SM) = d+{SM) U d-{SM), 

d±{SM) = {{x, v) G d{SM), T Wi^), ^) > }. 

The standard volume forms on SM and d{SM) that we will use are defined by 

ciE2"-i = dV'AdS^ 
dE2"-2 = dV'-^AdS^ 

where dV"' (resp. dV""'^) is the volume form of M (resp. dM ), and dSx = a/ det g{x)dE. 
where dEj. is the Euclidean volume form of S^ in T^M. For {x,v) G d{SM), let 
fi{x,v) = \{h'{x),v)\ and let L'j^{d^{SM)) be the space of functions on d^{SM) with 
inner product 

Jd+(SM) 

Without loss of generality, we may assume that (M, g) is embedded in (A^, g) where 
is a compact n-dimensional manifold without boundary. Let ipt be the geodesic flow 
on and X = -^Lp^lf^Q be the geodesic vector field. If {x,v) G SM, let 'y{t,x,v) be 
the unit speed A^-geodesic starting from x in the direction of v. Then 

Vt{x,v) = {-f{t,x,v),^{t,x,v)). 

Define the travel time r : SM — )• [0, oo] by 

t{x,v) = mf{t > : 7(t,a;,f) G N\M}. 

We say that (M, g) is non-trapping if r(x, f ) < oo for all {x, v) G SM. 

Definition. The geodesic ray transform of a function f G C°^{SM) is the function 

t(x,v) 

If{x,v)= j f{^tix,v))dt, {x,v)ed+{SM). 



Note that if the manifold (M, g) is non-trapping and has strictly convex boundary, 
then / : C°°(S'M) C{d+{SM)), and Santalo's formula [TO] implies that I is also a 
bounded map L'^{SM) — )■ L^(9+(S'M)). The general problem in tensor tomography is 
to determine properties of a function / from its integrals over geodesies as encoded by 
the transform If. 

Question. Given f G C°°{SM), what properties of f may be determined from the 
knowledge of I f ? 



TENSOR TOMOGRAPHY: PROGRESS AND CHALLENGES 



3 



Clearly a general function / on SM is not determined by its geodesic ray transform 
alone, since / depends on more variables than //. In applications one often encounters 
the transform I acting on special functions on SM that arise from symmetric tensor 
fields, and we will now consider this case. 

Let / = fi^...i^dx^'^ ® ■ ■ ■ ® (ix*" be a smooth symmetric m-tensor field on M. Such 
a tensor field induces a smooth function fm{x,v) on SM by 

The operator Im, defined by 

-^mf I fmi 

is called the geodesic ray transform of the symmetric tensor field /. If the manifold 
(M, g) is non-trapping and the boundary dM is strictly convex, then 

J„ : C^{M,S^{M)) ^ C{d+{SM)), 

where Sm{M) denotes the bundle of symmetric m-tensor fields over {M,g). We will 
frequently identify the tensor field f on M with the function fm on SM (see [32] for 
more details). 

It is known that any symmetric smooth enough tensor field / may be decomposed 
in a potential and solenoidal part : 

f = f + dp, sr = o, p\dM = o, 

where p is a smooth sjTiimetric (m — l)-tensor field on M, the inner derivative d = aV 
is the symmetric part of the covariant derivative V, and 6 is the divergence (the adjoint 
of —d in the natural inner product). If / is a 1-tensor, identified with a vector field 
W, this generalizes the usual Helmholtz decomposition of a vector field, 

W = W' + grad(p), div(Vr^) = 0, pIqm = 0. 

It is easy to see, using the fact that p vanishes on dM, that the geodesic ray transform 
of the potential part dp is zero. We denote by C^i^M, Sm{M)) the space of smooth 
solenoidal m-symmetric tensor fields. The remark above means that we can only expect 
to recover the solenoidal part of a tensor field from its ray transform. This leads to the 
following definition of solenoidal injectivity, or s-injectivity for short. 

Definition. The ray transform on symmetric m-tensors, m > 1, is said to be s- 
injective if Imf = implies /'^ = for any f G C°°{M, Sm{M)) . In the case of 
functions on M (m = 0), Iq is said to be s-injective if lof = implies / = for any 
f e C°^{M). 

The transforms Im arise in several applications as well as in the boundary rigidity 
problem. The latter consists in determining the Riemannian metric of a compact 
Riemannian manifold with boundary, modulo isometries fixing the boundary, from the 
distance function dg\gMxdM between boundary points [2S]- The case of /q when the 
metric is Euclidean is the standard X-ray transform that integrates a function along 
lines. Radon found in 1917 an inversion formula to determine a function knowing the 
X-ray transform. Inversion formulas of this type have been implemented numerically 
using the filtered backprojection algorithm which is used today in CT scans. 
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Another important transform in medical imaging and other apphcations is the Doppler 
transform which integrates a vector field along lines. This corresponds to the case of 
Ji for the case of the Euclidean metric. The motivation is ultrasound Doppler tomog- 
raphy. It is known that blood flow is irregular and faster around tumor tissue than 
in normal tissue and Doppler tomography attempts to reconstruct the blood flow pat- 
tern. Mathematically the problem is to what extent a vector field is determined from 
its integral along lines. 

The case of integration along more general geodesies arises in geophysical imaging in 
determining the inner structure of the Earth since the speed of elastic waves generally 
increases with depth, thus curving the rays back to the Earth surface. It also arises 
in ultrasound imaging. The geodesic ray transform Jq, that is, the integration of a 
function along geodesies, arises as the linearization of the boundary rigidity problem 
in a conformal class of metrics. The linearization of the boundary rigidity problem 
itself leads to I2, i.e. the integration of tensors of order two along geodesies. The case 
of integration of tensors of order 4 along geodesies arises in certain inverse problems in 
elasticity [42j . 

Many of the results in this survey are valid in the case when (M, g) is simple, a 
notion that naturally arises in the context of the boundary rigidity problem |26j. We 
recall that a Riemannian manifold with boundary is said to be simple if the boundary 
is strictly convex and if any two points are connected by a unique geodesic depending 
smoothly on the endpoints. In particular, a simple manifold is nontrapping and has 
no conjugate points. 

One of the main results we review in this paper is the s-injectivity of Im for all m 
for simple two-dimensional manifolds that was proven recently in [32j. 

Theorem 1.1. If {M,g) is a simple two-dimensional manifold, then Im is s-injective 
for any m > 0. 

This result was known earlier for m = [27] , m = 1 [3] and m = 2 ^6] . A key point 
in proving the result for general m is the efficient use of surjectivity properties of Jq, 
the adjoint of Jq. In fact, ^32] gave the following more general result. 

Theorem 1.2. // {M,g) is a compact non-trapping two dimensional manifold with 
strictly convex boundary, and if Iq and Ii are s-injective and Iq is surjective, then Im 
is s-injective for any m > 0. 

To describe in detail the adjoint Jq, for any function w on 9+(S'M) we define the 
function 

W^{X,V) = W{ip_rix-v){x,v)), {X,V) G SM. 

Then the solution of the boundary value problem for the transport equation 

Xu = in SM, u\g^(^sM) = w 

is equal to u = w^p. 

Recall that / is a bounded map L'^{SM) L'^^{d+{SM)). The adjoint J* is bounded 
L^(5_|_(5'M)) — )■ L'^{SM), and it is easy to compute explicitly. In the case of Jq, for 
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/ e C°°(M) and w G C°°{d+{SM)), we have 

p i't{x,v) 

iIof,w)Li(d+(SM)) = I / f{(ptix,v))w^{iptix,v))fxdtdJ:^''~^ 

d+(SM) Jo 
SM 

fix) ( / w4x,v)dS,{v)) 

The second equahty used Santalo's formula (TU]. From this computation we conclude 
that 

Iqw{x) = / w-^{x,v) dSxiv). 

J Sx 

Similarly, the adjoint of is the operator I^^ : {d+{SM)) (M, Sm{M)) which 

is given by 



Sx 

Definition. We say that Jg is surjective if for any f G C°^(M), there is a function 
w G C°^{d+{SM)) with I*w = f inM and G C°°(5M). 

The surjectivity of Jq in the above sense was proved in [32] on simple manifolds of 
any dimension. We will show below how this result is used in the uniqueness proof of 
tensor tomography in two dimensions. 

In this paper, we also review results in higher dimensions. Here is a summary of 
what is known about s-injectivity on simple manifolds of dimension n > 2: 

• Iq is injective [27] . 

• Ii is s-injective [3]. 

• Im is s- injective for all m if n = 2 [32]. 

• Im is s-injective for all m for manifolds of negative sectional curvature [37], or 
under certain other curvature restrictions [9], [36], [12] . 

• I2 is s-injective for generic simple metrics including real-analytic ones [52] . 

See [9], [13], [15], [53], [56] for uniqueness results on certain non-simple manifolds. We 
will also review results on the stability and range for J^, and moreover we propose 
several open problems. 

A brief summary of the contents of this paper is as follows. Section [2] contains 
preliminaries and notation used in the paper. In Section [3] we review the two proofs 
of Theorem 11.11 given in [32]. In Section H] we explain a natural approach to the proof 
of the so-called Pestov identities used in Section [31 This energy estimate approach 
resembles Carleman estimates. In Section [5] we review a microlocal approach to the 
study of the geodesic ray transform that gives in particular stability estimates which 
are summarized in Section In Section [7] we consider the scattering relation which 
is used in the characterization of the range and is of independent interest. In Section 
M we state the result of [35j on the range of the geodesic ray transform. In Section 
we summarize several results for the attenuated ray transform for unitary connections 
proved in [33]. In Section [10] we survey the result of [51] on s-injectivity of the ray 
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transform on 2-tensors on closed Anosov surfaces. Finally in Section [TT] we state several 
open problems. 

Acknowledgements. M.S. was supported in part by the Academy of Finland and 
an ERC starting grant, and G.U. was partly supported by NSF and a Walker Family 
Endowed Professorship. 

2. Facts about the unit circle bundle 

This section contains some facts needed for explaining the uniqueness proof for ten- 
sor tomography on surfaces, and we will restrict our attention to two dimensional 
manifolds. Let (M, g) be a compact oriented two dimensional Riemannian manifold 
with smooth boundary dM. As usual SM will denote the unit circle bundle which is 
a compact 3-manifold with boundary given by d{SM) = {{x,v) G SM : x G DM}. 

Let X denote the vector field associated with the geodesic flow iff Since M is 
assumed oriented there is a circle action on the fibers of SM with infinitesimal gen- 
erator V called the vertical vector field. It is possible to complete the pair X, to a 
global frame of T{SM) by considering the vector field X± defined as the commutator 
X_|_ := [X, V"]. There are two additional structure equations given by X = [V,X_|_] and 
[X, X_l] = —KV where K is the Gaussian curvature of the surface. Using this frame we 
can define a Riemannian metric on SM by declaring {X, X^, V} to be an ortho normal 
basis. This metric coincides with the Sasaki metric on SM, and the volume form of 
this metric will be denoted by rfS'^. The fact that {X, X^, V} are orthonormal together 
with the commutator formulas implies that the Lie derivative of dS^ along the three 
vector fields vanishes, in other words, the three vector fields preserve the volume form 
(iS^. See [?5] for more details on these facts. 

It will be useful to have explicit forms of the three vector fields in local coordinates. 
Since {M,g) is two dimensional, we can always choose isothermal coordinates (xi,X2) 
so that the metric can be written as ds"^ = e^^{dx\ + dx'^ where A is a smooth real- 
valued function of x = (xi,a;2). This gives coordinates (xi,X2,6') on SM where is 
the angle between a unit vector v and d/dxi. In these coordinates the vertical vector 
field is just 



and the other vector fields are given by 




, d . ^ d ( d\ . ^ 
+ sm 9— + -— sm 9 + 

OXi 0x2 \ 0x1 




X^ = -e-^ ~sm9— + cos9- — cos^ + — sm^ — 

\ 0x1 0x2 \oxi 0x2 J o9 

Given functions m, v : SM — ^ C we consider the inner product and norm 

uv dT? , \\u\\ = [u^uY^"^ . 

SM 
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Since X, V are volume preserving we have {Vu, v) = —{u, Vv) for u, v G C°°{SM), 
and if additionally u\q(^sm) = or v\q(^sm) = then also {Xu,v) = —{u,Xv) and 
{X±u,v) = -{u,Xj_v). 

The space L'^{SM) decomposes orthogonally as a direct sum 

L\SM) = ^H, 

fcez 

where is the eigenspace of —iV corresponding to the eigenvalue k. A function 
u G L'^{SM) has a Fourier series expansion 

oo 

U= ^ Mfc, 
fc=— oo 

where Uk G Hk- Also = Yl ll'^fclP? where = {u,uY^'^. The even and odd parts 
of u with respect to velocity are given by 

u+ ■= ^ Uk, u_ := ^ Mfc. 

fc even k odd 

In the (x, ^^)-coordinates previously introduced we may write 

Uk{x, 6)= Q- u{x, t)e~'^' d?j e'^^ = Uk{x)e'^''. 

Observe that for k > 0, Uk may be identified with a section of the k-th tensor power of 
the canonical line bundle; the identification takes Uk into Uke-^^{,dzY where z = Xi+ix2- 
The next definition introduces holomorphic and antiholomorphic functions with re- 
spect to the 6 variable. 

Definition. A function u : SM C is said to be holomorphic if Uk = for all k < 0. 
Similarly, u is said to he antiholomorphic if Uk = for all k > 0. 

Let flk '■= Hk n C°°{SM). As in [12] we introduce the following first order elliptic 
operators 

given by 

77+ := (X + iX^)/2, r/_ := (X - iX^)/2. 

Clearly X = rj^ + ri_. From the structure equations for the frame {X, X±,V} one 
easily derives: 

1]+ -Mk ^ ^k+i, ri^ -Mk ^k~i, iv+y = -V-- 

We will also employ the fiberwise Hilbert transform H : C°°{SM) — )■ C°°{SM), 
defined in terms of Fourier coefficients as 

Huk ■■= -isgn{k)uk. 

Here sgn(fc) is the sign of k, with the convention sgn(O) = 0. Thus, u is holomorphic 
iff (Id — iH)u = Uq and antiholomorphic iff (Id + iH)u = uq. 

The following commutator formula for the Hilbert transform and the geodesic vector 
field, proved in [52], has been a crucial component for many results reviewed in this 
paper. 
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Proposition 2.1. Let {M,g) be a two dimensional Riemannian manifold. For any 
smooth function u on SM we have the identity 

[H,X]u = X^u^ + {X^u), 

where 

Uq{x) = — / u{x,v)dSx 



271 . _ 

is the average value. 

Proof. It suffices to show that 

[Id + iH, X]u = iX^uo + i{X^u)o. 

Since X = 77+ + ?7_ we need to compute [Id + iiY, ?7±], so let us find [Id + iH,ri+]u, 
where u = J2k '^k- Recall that (Id + iH)u = uo + 2 X]fc>i ""fc- '^^d: 

(Id + iH)r]+u = ?7+M_i + 2 ^ r]+Uk, 



k>0 



7]+ (Id + iH)u = 7]+uq + 2 ^ ?7+Mfc. 



k>l 

Thus 

[Id + iH, r]^]u = Tj+u-i + r]+Uo. 

Similarly we find 

[Id + iH, ri^]u = —77- Mo — ri^ui. 
Therefore using that iX± = rj^ — t]^ we obtain 

[Id + iH, X]u = iX±uo + i{X±u)o 

as desired. □ 

3. Tensor tomography on surfaces 

The paper [32] gave two proofs for uniqueness in tensor tomography on a simple 
surface (M, g). In this section we will give an outline of both proofs. They are based on 
Pestov identities, which are energy estimates for operators related to the ray transform, 
and which will be discussed in more detail in Section SI Below we will make use of the 
concepts introduced in Sections [T] and [2J 

First proof. To explain the idea behind the first proof of s-injectivity, let us first 
assume that / is a 0-tensor, that is, / G C°°{M). Assuming that Jq/ = 0, it is required 
to show that / = 0. The first step is a reduction from the integral operator Jq into a 
PDE question involving a transport equation. The function 

/•t{x,v) 

u{x,v) = I f{(pt{x,v))dt, {x,v) G SM 



solves the transport equation 

Xu = -/ in SM, u\o(^sM) = 0. 
It is enough to show that u = 0, since then also / = 0. 
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Isothermal coordinates allow to identify 

SM = {{x,e); xeB,ee [0,27r)}. 
The vertical vector field on SM is V = We want to show that 

I ■' ^ u = 0. 

U\d{SM) = U 

If / is a 0-tensor, / = f{x), then Vf = 0. Thus it is enough to show that 

VXu = 



I „ u = 0. 

u\a{SM) = U 

This calls for a uniqueness result for the operator P = VX. In isothermal coordi- 
nates, this operator has the form 

where h{x,d) is a certain smooth function. It turns out that the operator P is rather 
exotic and there do not seem to be general results on uniqueness properties of such 
operators in the literature. Here are some facts about the operator P: 

• it is a second order operator on 3D manifold SM 

• it has multiple characteristics 

• P + W has compactly supported solutions for some first order perturbation W 

• it enjoys a subelliptic type estimate ||'u||i^i(5Af) < C||-Pw||i2(5jvj) for u G C°°{SM) 
with u\d(^sM) = 0. 

However, we can still prove a global uniqueness result for P by using energy esti- 
mates. This involves the Pestov identity in L'^{SM) inner product when u\q(^sm) = 0: 

\\Puf = \\Auf + \\Buf + {i[A,B]u,u) 

where P = A + iB, A* = A, B* = B. 

We will compute the commutator below, and this gives (see Proposition 14. 2p 

\\Pu\\^ = \\XVu\\^ - {KVu, Vu) + \\Xuf. 

It is known [53] that on simple manifolds 

\\XVuf - {KVu, Vu) > 0, ue C^{SM),u\9^sM) = 0. 

(Note that in the case of non-positive curvature, i.e. K < 0, one always has ||XV^Mp — 
{KVu, Vu) > 0.) Thus Pu = implies u = 0, showing injectivity of Jq. 

We now return to tensor tomography. Let Xu = —f in SM, u\q(^sm) = where 
/ is the function on SM corresponding to a symmetric m-tensor field. It will be 
convenient to switch to a slightly different setup and think of u and / (which are 
functions SM — )■ C) as sections of the trivial bundle E = SM x C. The transport 
equation then becomes an equation for sections of E, 

D\u = -f 

where Dx = d is the fiat connection on the trivial bundle E. 
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One benefit of tfiis (trivial) change of point of view is that from the equation on 
sections, one sees that the transport equation has a natural gauge group acting via 
multiplication by smooth functions c G C°°(M). This action preserves m-tensors, and 
leads to gauge equivalent equations 

Di{cu) = -cf 

where = d + A is a. gauge equivalent connection on E and A = —c^^dc is the 1-form 
determining the connection. 

Now we try to use an energy identity for the connections D"^. This Pestov iden- 
tity with a connection is proved in the same way as the usual Pestov identity (see 
Proposition 14. 3p . and reads in LF'{SM) norms 

\\V{X + A)Mf = ||(X + A)Vu\\'^ - {KVu, Vu) + ||(X + A)uf + {*FaVu,u). 

Here * is Hodge star and 

FA = dA + AAA 

is the curvature of the connection = d + A. We observe that if the curvature *Fa 
and the expression {Vu, u) have suitable signs, we gain a positive term in the energy 
estimate. 

This observation does not immediately lead to anything new since curvature is pre- 
served under gauge transformation. Thus, if is gauge equivalent to , then 
Fa = Fq = 0. However, we can use a generalized gauge transformation that arranges 
a sign for Fa- This involves gauge transformations via functions c that may depend 
on the V variable. Such transformations break the m-tensor structure of the equation, 
but turn out to be manageable if the gauge transforms are holomorphic in a suitable 
sense. 

Recall from Section [2] that a function u G L'^{SM) is called holomorphic if Uk = 
for k < 0. The main point is the following theorem guaranteeing that holomorphic 
gauge transformations always exist. This is related to injectivity of the attenuated ray 
transform on simple surfaces [41j, and in the form below it is proved in [32j and [33j. 
The proof is based on the surjectivity of Iq. 

Theorem 3.1 (Holomorphic gauge transformation). If A is a 1-form on a simple 
surface, there is a holomorphic w G C°^(S'M) such that X-\-A = e'^oXo . 

Proof. Since M is simply connected, there is a Hodge decomposition Aj dx^ = da + -kdh 
for some a,b & C°°{M) (★ is the Hodge star operator). In terms of the corresponding 
functions on SM we have A = Xa + Xj_b. Replacing w hj w — a, it is enough to 
consider the case where A = Xj_b. 

Let us try a solution of the form w = (Id + iH)w where w G C°°{SM) is even with 
respect to v. By Proposition 12. ![ 

Xw = (Id + iH)Xw - i[H, X]w = (Id + tH)Xw - iX±wo. 

Now it is sufficient to find w even with Xw = and Wq = —ib. Using the surjectivity 
of /q I39j, there is some h G C~(a+(5M)) with I*h = -2nib. But if w' G C^{SM) is 
the function with Xw' = in SM and w'\q^(^sm) = we have (w')o = 2^-^o^ = 
It is enough to take w to be the even part of w' with respect to f . □ 
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We can now explain the end of the proof of the uniqueness result for tensor tomog- 
raphy on simple surfaces. Let / = XlfcL-m fk be an m-tensor written in terms of its 
Fourier components, and let 

Xu = -/, u\d(^sM) = 0. 

Choose a primitive if of the volume form Ug of {M,g), so that dip = Ug. Let s > be 
large, let Ag = —isip, and choose a holomorphic w with X + As = e"^^ o X o e"^^ . The 
transport equation becomes 

(X + As){e''-u) = -e'^f, e'^u\9isM) = 0. 

Here the curvature of Ag has a sign and one has information on Fourier coefficients 
of e*""/. The Pestov identity with connection allows to control Fourier coefficients of 
e'^^u, eventually proving s-injectivity of Im- 

Heuristically, the proof above involves "twisting" the trivial bundle E hj a holo- 
morphic gauge transformation to make it positively curved, using the Pestov identity 
with a large positive term coming from the connection to absorb error terms, and then 
undoing the gauge transformation (this is possible because of holomorphicity) to get 
uniqueness. This idea of twisting to impose positivity to prove a vanishing theorem is 
of course well known in Complex Geometry and it is the way one proves results like 
the Kodaira vanishing theorem [171 . Our setting is more complicated since the relevant 
PDE is the transport equation which is harder to handle than the Cauchy-Riemann 
equation. However this analogy is important and permeates all work; in particular 
the injectivity results on the attenuated ray transform for unitary connections, to be 
discussed later on, are also proved in this fashion. 

There is an interesting connection between the Pestov identity with connection Ag 
above and with Carleman estimates. In fact, the Pestov identity with Ag implies the 
estimate 

Here we use the norms 




Formally this looks very much like a Carleman estimate with exponential weights, but 
it involves some slightly exotic spaces and one can see that the positivity comes from 
lm{w) (not Re(w) as is usual in Carleman estimates)! We finally remark that such an 
estimate is sufficient for 

• absorbing large attenuation (even for systems, see Section ^ 

• absorbing error terms coming from m-tensors. 

However, it seems that the estimate may not be enough to 

• localize in space 

• absorb error terms coming from curvature of M. 

Second proof. Next we explain a very short alternative proof to a key step in the 
injectivity result. 
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Suppose that m is a smooth solution of Xu = — / in SM where fk = for k < —m — l 
and u\q(^sm) = 0. We wish to show that = for /c < — m. This, together with the 
analogous result for positive Fourier coefficients, implies that / = Xh where the Fourier 
expansion of h has degree m — l and h\Q(^sM) = 0, thus proving s-injectivity. 

We choose a nonvanishing function h G Qm- In fact, in isothermal coordinates, we 
can set 

h{x,y,e):=e'"^'. 

Define the 1-form 

A := -h-^Xh. 

Then hu solves the problem 

{X + A){hu) = -hf in SM, hu\9^sR4) = 0. 

Note that hf is a holomorphic function. Next we employ a holomorphic integrating 
factor, as above: by Theorem 13.11 there exists a holomorphic w G C°°{SM) with 
Xw = A. The function e^hu then satisfies 

Xie'^hu) = -e'"hf in SM, e'"hu\d(^sM) = 0. 

The right hand side e^hf is holomorphic. It is known that the solution e^hu, which 
vanishes on d{SM) also has to be holomorphic and further [e^huj^ = 0. This follows 
from the s-injectivity of /q and Ji (see [32], [S]). Looking at Fourier coefficients shows 
that {hu)k = for < 0, and therefore Uk = for k < —m as required. 

4. Pestov identity 

In this section we consider the Pestov identity, which is the basic energy identity 
that has been used since the work of Mukhometov [27] in most injectivity proofs of 
ray transforms in the absence of real-analyticity or special symmetries. Pestov type 
identities were also used in [3] to prove s-injectivity for Ji on simple manifolds and in 
[2Z] to prove s-injectivity for any m in any dimensions if the sectional curvatures are 
negative. See [H], [SS], [12] for further results. Pestov identities have often appeared 
in a somewhat ad hoc way, but here we follow [32] which gives a new point of view 
making the derivation of these identities more transparent. We will only consider two 
dimensional manifolds in this section. 

The easiest way to motivate the Pestov identity is to consider the injectivity of the 
ray transform on functions. The first step, as discussed in Section ??, is to recast the 
injectivity problem as a uniqueness question for the partial differential operator P on 
SM where 

P := VX. 

This involves a standard reduction to the transport equation. 

Proposition 4.1. Let {M,g) be a compact oriented nontrapping surface with strictly 
convex smooth boundary. The following statements are equivalent. 

(a) The ray transform I : C°°{M) — )■ C{d^{SM)) is injective. 

(b) Any smooth solution of Pu = in SM with u\d(^sM) = is identically zero. 
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Proof. Assume that the ray transform is injective, and let u G C°°{SM) solve Pu = 
in SM with u\d(SM) = 0. This implies that Xu = —f in SM for some smooth / only 
depending on x, and we have = M|a+(SM) = If- Since / is injective one has / = 
and thus Xu = 0, which implies m = by the boundary condition. 

Conversely, assume that the only smooth solution of Pu = in SM which vanishes 
on d{SM) is zero. Let / G C°^(M) be a function with If = 0, and define the function 

pt{x,v) 

u{x,v) := / /(7(t, X, f )) dt, {x,v) E SM. 
Jo 

This function satisfies the transport equation Xu = — / in SM and u\q(^sm) = since 
// = 0, and also u G C°°{SM) (see [33j). Since / only depends on x we have Vf = 0, 
and consequently Pu = in SM and u\q(^sm) = 0. It follows that u = and also 
/ = -Xu = 0. □ 

We now focus on proving a uniqueness statement for solutions of Pu = in SM. 
For this it is convenient to express P in terms of its self-adjoint and skew-adjoint parts 
in the L'^{SM) inner product as 

p _|_ TD* TD p* 

P = A + iB, A:=—- , B := . 

Here the formal adjoint P* of P is given by 

P* := XV. 

In fact, if M G C°°{SM) with u\a(sM) = 0, then 

(1) \\Puf = {{A + tB)u, {A + iB)u) = \\Auf + \\Buf + i{Bu, Au) - t{Au, Bu) 

= \\Auf + \\Buf + {i[A, B]u, u). 

This computation suggests to study the commutator i[y4,i?]. We note that the argu- 
ment just presented is typical in the proof of Carleman estimates [21j. 

By the definition of A and B it easily follows that 5] = i[P*,P]. By the 
commutation formulas for X, X_^_ and V , this commutator may be expressed as 

[P\ P] = xvvx - vxxv = vxvx + x^vx - vxvx - vxx^ 

= V[X^, X]~X^ = -X^ + VKV. 

Consequently 

{[P*, P]u, u) = \\Xuf - [KVu, Vu). 

If the curvature K is nonpositive, then [P*, P] is positive semidefinite. More generally, 
one can try to use the other positive terms in ([1]). Note that 

\\Au\\^ + \\Buf = ^{\\Pu\\^ + \\P*uf). 

The identity ([1]) may then be expressed as 

\\Puf = \\P*uf + {[P*,P]u,u). 

(Note that we could have just started from the last identity, but expressing matters via 
A and B highlights the similarity to Carleman estimates.) Moving the term H-PmH^ to 
the other side, we have proved the version of the Pestov identity which is most suited 
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for our purposes. The main point in this proof was that the Pestov identity boils down 
to a standard estimate based on separating the self-adjoint and skew-adjoint parts 
of P and on computing one commutator, [P*, P]. 

Proposition 4.2. If{M,g) is a compact oriented surface with smooth boundary, then 

\\XVuf - {KVu, Vu) + \\Xuf - \\VXuf = 

for any u G C°°{SM) with u|a(5Af) = 0. 

It is known [13], [23] that on a simple surface, one has 

\\XVuf - {KVu, Vu) > 0, ue C°°{SM), u\d(^sM) = 0. 

Also, if Xu = — / where / = /o + /i + f~i is the sum of a 0-form and 1-form, we have 

\\Xuf - \\VXuf = ll/of > 0. 

These two facts together with the Pestov identity give the standard proof of s-injectivity 
of the ray transform for 0-forms and 1-forms on simple surfaces. It is easy to see where 
this proof breaks down if m > 2: the Fourier expansion / = YlT=-m fk implies 

\\xur-\\vxur = \\for- ^ (e-mf.r. 

2<\k\<m 

This term may be negative, and the Pestov identity may not give useful information 
unless there is some extra positivity like a curvature bound. 

Finally, we consider the Pestov identity in the presence of attenuation given by 
A{x,v) = Aj{x)v^ where Aj dx^ is a purely imaginary 1-form on M. We write A both 
for the 1-form and the function on SM . The geometric interpretation is that d + Ais 
a unitary connection on the trivial bundle M x C, and its curvature is the 2-form 

Fa:= dA + A ^ A. 

Then -kF^ is a function on M where -k is the Hodge star. We consider the operator 

p ■= V{X + A). 

Since A = —A, the formal adjoint of P in the L'^{SM) inner product is 

P* = {X + A)V. 

The same argument leading to Proposition 14. 2[ based on computing the commutator 
[P*,P], gives the following Pestov identity proved also in [S3T, Lemma 6.1]. 

Proposition 4.3. // {M,g) is a compact oriented surface with smooth boundary and 
if A is a purely imaginary 1-form on M, then 

\\{X + A)Vu\\'^ - {KVu,Vu) + ||(X + A)Mf - \\V{X + A)u\\'^ + {i.FAVu,u) = 

for any u G C°°(5'M) with u\q(sm) = 0. 

Using the Fourier expansion of m, the last term in the identity is given by 

oo 

ik{'kFAUk,Uk) 
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This shows that if u is holomorphic and z * > 0, or if u is antiholomorphic and 
i-k Fa < 0, one gains an additional positive term in the Pestov identity. This is crucial 
in absorbing negative contributions from the term + — || V^(X + 74)m|P when 

proving s-injectivity on tensor fields. 



5. MiCROLOCAL APPROACH 

A different approach that is useful to prove s-injectivity of Im in some cases and 
gives stability estimates as well as reconstruction formulas in some cases was started in 
[50] and developed further in [511 [52], [531 Ell- We describe the method in more detail 
for /q. Let {M,g) be a simple manifold embedded in a closed manifold {N,g) and let 
f/ be a simple neighborhood of M in N. 

Theorem 5.1. Jq/q is an elliptic pseudodifferential operator on U of order -1. 
Proof. It is easy to see, that 

t(x,v) t{x,v) 

/(7 {t,x,v))dt. 

Sx —t{x, — v) Sx 



(2) (/*/o/) (x) = jdS, j / (7 (t, x, v)) dt = 2jdS, J 



Before we continue we make a remark concerning notation. We have used up to 
now the notation 'j{t,x,v) for a geodesic. But it is known, that a geodesic depends 
smoothly on the point x and vector ^t G Tj,.(M). Therefore in what follows we will also 
use sometimes the notation 7(x, vt) for a geodesic. Since the manifold M is simple any 
small enough neighborhood U (in (A^, g)) is also simple (an open domain is simple if its 
closure is simple). For any point x E U there is an open domain C [U) such 
that exponential map expx '■ — )■ U, exp^rj = 7(x, rj) is a diffeomorphism onto U. Let 
Dx, X G M be the inverse image of M, then expx{Dx) = M and expx\Dx '■ Dx — ^ M 
is a diffeomorphism. 

Now we change variables in (2), y = 7(x, vt). Then t = dg {x, y) and 

(/*//) (x) = J K{x,y) f{y)dy, 

M 

where 

det (exp-^)' (x, y) ^/detgJx) 



K{x,y) = 2- 



dg ^{x,y) 
Notice, that since 

(3) -f{x,ri) = X + r] + 0{\rif), 

it follows, that the Jacobian of the exponential map is 1 at 0, and then det{exp~^ )'(x, x) 
1/ det (expx)' {x, 0) = 1. From (3) we also conclude that 

d^ (x, y) = G^j (x, y) (x - yf (x - yY , (x, x) = g^j (x) , G,^ e {M x M) . 
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Therefore the kernel of Jq/q can be written in the form 

2 det {exp~^)' (x, y) \/ det g (x) 



Gij (x, y) {x - yf {x - yf^ 



{n-l)/2- 



Thus the kernel K has at the diagonal x = y a singularity of type |a; — y| . The 
kernel 

I N 2^detc/(s) 

^0 y) = ^ (n-l)/2 

9ij {x) {x - y)* (x - y)^ 



has the same singularity. Clearly, the difference K — Kq has a singularity of type 
\x — • Therefore the principal symbols of both operators coincide. The principal 

symbol of the integral operator, corresponding to the kernel Kq coincide with its full 
symbol and is easily calculated. As a result 

f e"*(^''^) _i 

(t(/*/o) = 2Vdetc/(x) / -——-T—[^Hdy = Cn\i\ ■ 

J {9ij{x)y'y^y 

□ 

Let g he a. simple metric in M. Extend g near M and let Mi be a simple manifold 
with boundary so that M is a compact subset of Mi. We will work with / supported 
in M. We assume that / is extended as outside M. Choose a smooth function x 
supported in Mi such that x = 1 iisar M. 

It was shown in [51] that the normal operator A^^^ = I^Im is a pseudodifferential 
operator of order -1, for m = 0,1,2 which is elliptic acting on solenoidal tensor fields. 
We have. 

Theorem 5.2. Let g be a simple metric in M and let x be as before. Then one can 
construct a pseudodifferential operator aijki{x, D) of order 1 so that for any symmetric 
2-tensor f G L'^{M) we have, 

(4) xa^M^, D)x{Njf' = fl,^ + Kf, 

where K : L'^{M) — )■ if^(Mi) is bounded. Here /|^^ denotes the solenoidal part of f on 
Ml. 

This result was extended to tensor fields of any order in [17]. 

When (7 is a real-analytic simple metric it was shown in [52] that I2 is s-injective. The 
proof constructs a parametrix as in the previous result with K analytic regularizing, 
that is Kf is real-analytic on Mi for / G L'^{M). The idea of the proof for /q is that 
if Jo/ = 0, / G L'^{M) then / = —Kf. Since Kf is real-analytic on Mi and supported 
on M it must be zero. For the details of the proof for I2 see [52]. 

6. Stability estimates 

It was shown in \5T], [17] that for a simple manifold s-injectivity of Im implies stability 
estimates. This is based on the fact that Ng = I^Im is an elliptic pseudodifferential 
operator acting on solenoidal tensor fields. We have the following stability estimate for 

/o my- 
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Theorem 6.1. Let g be a simple metric in M and assume that g is extended smoothly 
as a simple metric near the simple manifold Mi DD M. Then for any function f G 



|l2(a/)/C < \\Ngf\\Hi(Mi) < C||/||l2(m). 

Similarly s-injectivity of Ji implies the stability estimate: 

Theorem 6.2. Assume that g is simple metric in M and extend g as a simple metric 
in Ml DD M. Then for any 1-form f = fidx^ in LP'{M) we have 

A sharp stability estimate for I2, assuming that I2 is known to be s-injective, was 
proved in 



Theorem 6.3. Let g he a simple metric in M and assume that g is extended smoothly 
as a simple metric near the simple manifold Mi DD M . Also assume that I2 is s- 
injective. Then for any symmetric 2-tensor field f in L'^{M), 

\\f'\\L^(M)/C < ||iVg/||/fi(Mi) < C'||/'*||l2{M)- 

In order to describe possible stability estimates for we describe an earlier result 
for I2. In order to state the result we first take boundary normal coordinates x^, ...,x^ 
with x"' = the defining function of dM. Introduce the space H^{Mi) with norm equal 
to the norm outside a neighborhood of dM and near dM (but outside M) having 
the following form in normal local coordinates: 

(5) i/iii(Mo=/ (Ei^^/r+K^n/r+i/r supp/c^/. 

"'A^i \i=i J 

Here [/ is a small neighborhood of a point on dM and the norm in H^{Mi) is defined 
by using a partition of unity. 
Next we define the norm 

n 
i=l 

The earlier stability result for I2 is: 

Theorem 6.4. Assume that g is simple metric in M and extend g as a simple metric 
m i Ml DD M. 

(a) The following estimate holds for each symmetric 2-tensor f in H^[M): 

liriL.(M) < c||a^p/IIh2(mo + Csll/llH-^(AfO' V. > 0. 

(h) Ker/2 n SL'^{M) is finite dimensional and included in C°°{M). (S stands for 
solenoidal). 

(c) Assume that I2 is s-injective in M , i.e., that Ker/^ fl SL'^{M) = {0}. Then for 
any symmetric 2-tensor f in H^{M) we have 

||/''|lL2(jvf) < C\\Ngf\\jj2(^f,f^y 
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This result was proven in [51] for the case m = 2. Using the resuhs of [47], stabihty 
estimates of this type can be shown to be vahd for any m. 

7. The scattering relation 

To state the results for the range of for simple surfaces we need to recall the 
definition of the scattering relation which is a subject of interest in its own right. 

Suppose we have a Riemannian metric in Euclidean space which is the Euclidean 
metric outside a compact set. The inverse scattering problem for metrics is to deter- 
mine the Riemannian metric by measuring the scattering operator (see [iSj). A similar 
obstruction to the boundary rigidity problem occurs in this case with the diffeomor- 
phism ip equal to the identity outside a compact set. It was proven in [TH] that from the 
wave front set of the scattering operator, one can determine, under some non-trapping 
assumptions on the metric, the scattering relation on the boundary of a large ball. 
This uses high frequency information of the scattering operator. In the semiclassical 
setting Alexandrova has shown for a large class of operators that the scattering op- 
erator associated to potential and metric perturbations of the Euclidean Laplacian is 
a semiclassical Fourier integral operator quantized by the scattering relation [2]. The 
scattering relation maps the point and direction of a geodesic entering the manifold to 
the point and direction of exit of the geodesic. 

We proceed to define in more detail the scattering relation. To do this, let r° = 
Ad{SM) and note that this function is equal to zero on d^{SM) and is smooth on 
d+{SM). Its odd part with respect to f , 

t'_{x,v) = \{t\x,v)-t\x,-v)), 
is a smooth function on d{SM) (see for instance [TO]). 

Definition. Let (M, g) be non-trapping with strictly convex boundary. The scattering 
relation a : d{SM) — )■ d{SM) is defined by 

a{x,v) = (7(x, f , 2r° (x, f )), 7(x, f , 2r° (x, f ))). 

The scattering relation is a diffeomorphism d{SM) — )■ d{SM). Notice that a\d+{sM) '■ 
d+{SM) d-{SM), a\a_sM : d^{SM) ^ d+{SM) are diffeomorphisms as well. The 
manifold of inner vectors 9-|_(S'M) and outer vectors d-{SM) intersect at the set of 
tangent vectors 

do{SM) = {{x,v) G d{SM), {v{x),v) = }. 

Obviously, a is an involution, = id and do{SM) is the hypersurface of its fixed 
points, a{x,v) = (x,v), {x,v) G do^SM). 

A natural inverse problem is whether the scattering relation determines the metric 
g up to an isometry which is the identity on the boundary. This information takes 
into account all the travel times, not just the first arrivals like the boundary distance 
function. 

We remark that in the case that (M, g) is a simple manifold, and we know the metric 
at the boundary (and this is determined if dg is known), knowing the scattering relation 
is equivalent to knowing the boundary distance function ([26j). 
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We introduce the operators of even and odd continuation with respect to a: 

A±w{x,v) = w{x^v), (x, f ) G S+SM, 

Aj^w{x,v) = ± {a*w) (x,f), {x,v) e d^{SM). 
We will examine next the boundedness properties of A_, 

Lemma 7.1. A±: Ll{d+{SM)) -> Lf^^{d{SM)) are bounded. 
Proof. 

ld+(SM) Jd^{SM) 
d+{SM) Jd+{SM) 

where a : d+{SM) — )• d-{SM) is a diffeomorphism. Thus it is enough to show that 

Let w G C^id+iSM)). Then 

WTfidJ:^'"'^= / w^{^t{x,v))^idtdY?''-^ = w^^dT?""-^ 

d+{SM) Jd+(SM) Jo JsM 

Set u{x, v) = u{x, —v) for u G C°°{SM), one has 

->2n— 1 / „7, ^v'2n— 1 



a_(SM) Jo 

^ / 



2?i--2 



w;(a(?/,r/))(-/i)c/trfS2"-2 

Varying w shows that a*(-/idS2"-2) = /idS^"-^ 9+(5M)\ao^M. □ 
The adjoint : Lf^^{d{SM)) -> L2(9+(SM)) satisfies 

Jd+{SM) Jd^(SM) 
= I w(M±MOa)/ic/S2"~2 

ia+(5M) 

so A^U = (u ± M O a)\3^(^sM)- 

In |39j the following characterization of the space of smooth solutions of the transport 
equation was given. Here we define 

C^{d+{SM)) = {we C^{d+{SM)) : e C°°{SM)}. 

Lemma 7.2. 

C^{d+{SM)) = {we C^{d+{SM)) : A+w g C^{d{SM))}. 
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Then I*w e C^{M) whenever w e C^{d+{SM)). 

We conclude this section by defining certain operators which combine the operators 
A± introduced above with the fibrewise Hilbert transform H. These operators will 
be essential to determine the range of the ray transform in the next section. Set 
H±u = Hu± where (resp. U-) denote the even (resp. odd) part of m G C°^(S'M). 

We define 

(6) P_ = A*_H^A+, P+ = A*^H+A+. 

8. Range of the geodesic ray transform 

We now give the characterization of the range of Jq and Ji in terms of the scat- 
tering relation only. We have that these are the projections of the operators P_,P+ 
respectively (defined in ([H])). For the details see [SB]. 

Theorem 8.1. Let {M,g) be simple two dimensional compact Riemannian manifold 
with boundary. Then 

(1) A function u G {d^{SM)) belongs to the range of Iq iff u = P^w where 
weC^{d+{SM)). 

(2) A function u G {dj^{SM)) belongs to the range of Ii iff u = P+w where 
w G C^{d+{SM)). 

We now move on to describe the range of the geodesic ray transform for tensors of 
order > 2. For this we apply the ideas of the second proof of Theorem 11.11 described 
in Section 3. For the details see [35] . 

Let (M, g) be a simple surface. The metric g induces a complex structure on M 
and let k, be the canonical line bundle (which we may identify with T*M). Recall that 
Hm ijn G Z) is the set of functions in / G L^(S'M, C) such that Vf = imf . The set 
VLm = Hmf^ C^iSM, C) can be identified with the set r(M, k®™) of smooth sections 
of m-th tensor power of the canonical line bundle k. This identification depends on the 
metric and is explained in detail in [Ml Section 2], but let us give a brief description 
of it. Given a section ^ G r(M, k®'") we can obtain a function on simply by 
restriction to SM: ^ determines the function SM 3 {x,v) t— )■ ^x{v'^'^) and this gives a 
1-1 correspondence. 

Since M is a disk, there is ^ G r(M, k) which is nowhere vanishing. Having picked 
this section we may define a function h : SM — )■ by setting h{x,v) = ^x{v)/\C,x{v)\. 
By construction h & fli. Our description of the range will be based on this choice of 
h. Define 

Air^g = A:= -h-^Xh. 

Observe that since h G ^i, then = /t G Vt^i. Also Xh = rj^h + r]_h & 0.2® 
which implies that A &VLi ©fi_i. It follows that A is the restriction to SM of a purely 
imaginary 1-form on M, hence we have a unitary connection (see Section [9]). 

First we will describe the range of the geodesic ray transform I restricted to Vim'- 

:= : -> C°°(9+(SM), C). 

Observe that if u solves the transport equation Xu = —f with u\q^(^sm) = 0, then h~''^u 
solves (X - mA){h-"'u) = -h-"^ f and /i-™M|a_(5A/) = 0. Also note that /i"'"/ G Oq. 
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Thus 

where the left hand side is the attenuated ray transform of the unitary connection 
—mA. Attenuated transforms will be described in more detail in the next section, but 
the upshot is that we can prove a theorem similar to Theorem 18.11 but introducing this 
time a unitary connection as attenuation. Putting everything together one obtains a 
description of the range for as follows. Let 

/ w{x,v) ii{x,v)ed+iSM) 
(^^w[x,v). I (e— /;*^"'^A(^*(->-))'^t«;)oa(x,t;) ii (x^v) E d^SM) 

and 

In other words: 

/ w{x,v) ii{x,v)ed+{SM) 
L^mW{x, V) I (e-m/i(A)y^^ ^ a{x, v) if (x, v) e d.{SM) 

and 

Bmg=[e"'''^''\goa)-g]\9^^SM). 

We define 

Pm,— • BrnH —Qm- 

Then: 

Theorem 8.2 Let (M, g) be a simple surface. Then a function u e C°°{d^{SM), C) 

belongs to the range of Im if and only ifu = (/i™'|a+(5M)) Pm-W forw E S^{d+{SM), C), 
where this last space denotes the set of all smooth w such that QmW is smooth. 

Suppose now F is a complex-valued symmetric tensor of order m and we denote its 
restriction to SM by /. Recall from [321 Section 2] that there is a 1-1 correspondence 
between complex-valued symmetric tensors of order m and functions in SM of the 
form / = 'Y^^=-m fk where fk E Vtf. and fk = for all k odd (resp. even) if m is even 
(resp. odd). 

Since 

m 

lif) = J2 Wk) 

k=~m 

we deduce directly from Theorem 18.21 the following. 

Theorem 8.3. Let {M,g) be a simple surface. If m = 21 is even, a function u E 
C°°{d+{SM),<C) belongs to the range of the ray transform acting on complex-valued 
symmetric m-tensors if and only if there are W2k £ '5^((9+(5'M), C) such that 



u 



^ {h^^\d+{SM)) P2k-W2k- 



22 



G.P. PATERNAIN, M. SALO, AND G. UHLMANN 



Similarly, if m = 21 + 1 is odd, a function u G C°°(9+(5'M), C) belongs to the range of 
the ray transform acting on complex-valued symmetric m-tensors if and only if there 
are W2k+i G '5^+i(5+(S'M), C) such that 

I 

U= ^ {h'^''^^\d+{SM)) P2k+l-'W2k+l- 
k=-l-l 

9. Attenuated ray transform for unitary connections 

In this section we describe in detail certain injectivity results for the attenuated ray 
transform of a unitary connection [33]. We saw the appearance of the attenuated ray 
transform in the last section when we discussed the range of the (unattenuated) ray 
transform on tensors of any order. We also saw how useful was for the tensor tomogra- 
phy problem to introduce a connection to gain positivity in the Pestov identity. Here 
we take a closer and more systematic look. We motivate this section by discussing 
first another natural inverse problem: determine a unitary connection from its scat- 
tering relation, that is, parallel transport along geodesies between boundary points. 
Our results are for simple surfaces, but the definitions can be given in the context of 
non-trapping manifolds (M, g) with strictly convex boundary. 

Suppose E ^ M is a. Hermitian vector bundle of rank n over M and V is a unitary 
connection on E. Associated with V there is the following additional piece of scattering 
data: given {x,v) G d+{SM), let P{x,v) = P^{x,v) : E{x) — >• E{t{ o a{x,v)) denote 
the parallel transport along the geodesic 7(t,x,f). This map is a linear isometry and 
the main inverse problem we wish to discuss here is the following: 

Question. Does P determine V? 

The first observation is that the problem has a natural gauge equivalence. Let be a 
gauge transformation, that is, a smooth section of the bundle of automorphisms AvXE. 
The set of all these sections naturally forms a group (known as the gauge group) which 
acts on the space of unitary connections by the rule 

(^*V)s := %l)V{il)-^s) 

where s is any smooth section of E. If in addition iI)\qm = Id, then it is a simple 
exercise to check that 

-Pv = P%p*v- 

Thus we can rephrase the question above more precisely as follows: 

Question I. Let Vi and V2 be two unitary connections with Pvi = -Pv2- Does there 
exist a gauge transformation -0 with ipldu = Id and V'*Vi = V2? 

It is easy to see from the definition that a simple manifold must be diffeomorphic to 
a ball in M". Therefore any bundle over such M is necessarily trivial and from now on 
we shall assume that E = M x C^. 

Question I arises naturally when considering the hyperbolic Dirichlet-to-Neumann 
map associated to the Schrodinger equation with a connection. It was shown in that 
when the metric is Euclidean, the scattering data for a connection can be determined 
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from the hyperbolic Dirichlet-to-Neumann map. A similar result holds true on simple 
Riemannian manifolds: a combination of the methods in [16j and |55j shows that the 
hyperbolic Dirichlet-to-Neumann map for a connection determines the scattering data 
Pv. 

Elementary background on connections. Consider the trivial bundle M x C". 
For us a connection A will be a complex nxn matrix whose entries are smooth 1-forms 
on M. Another way to think of A is to regard it as a smooth map A : TM — )■ C"^"" 
which is linear in f G T^M for each x G M. 

Very often in physics and geometry one considers unitary or Hermitian connections. 
This means that the range of A is restricted to skew- Hermitian matrices. In other 
words, if we denote by u(n) the Lie algebra of the unitary group U{n), we have a 
smooth map A : TM — )■ u{n) which is linear in the velocities. There is yet another 
equivalent way to phrase this. The connection A induces a covariant derivative on 
sections s G C°°{M, C") by setting (Ias = ds + As. Then A being Hermitian or unitary 
is equivalent to requiring compatibility with the standard Hermitian inner product of 
C" in the sense that 

d{si, S2) = {dASi, S2) + (Si, dAS2) 

for any pair of functions Si,S2- 

Given two unitary connections A and B we shall say that A and B are gauge equiv- 
alent if there exists a smooth map u : M U{n) such that 

(7) B = u'^du + u'^Au. 

It is easy to check that this definition coincides with the one given in the previous 
section if we set ip = u~^. 

The curvature of the connection is the 2-form Fa with values in u{n) given by 

FA:=dA + AA A. 

If A and B are related by ([7j) then: 

Fb = Fa u. 

Given a smooth curve 7 : [a, h] — M , the parallel transport along 7 is obtained by 
solving the linear differential equation in C^: 

s + A(7(t),7(t))s = 0, 
s{a) = w G C". 

The isometry Pa{,i) : C" — )■ C" is defined as Pa{i){w) := s{h). We may also consider 
the fundamental unitary matrix solution U : [a,b] — t- U{n) of (jH]). It solves 



(9) 

Clearly Pa{i){w) = U{b)w. 



f7 + A(7(t),7W)f/ = 
U{a) = Id. 



The transport equation and the attenuated ray transform. Consider now the 
case of a compact simple Riemannian manifold. We would like to pack the information 
provided by (Q along every geodesic into one PDE in SM. For this we consider the 
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vector field X associated with the geodesic flow 
Ua:SM^ U{n) of 



and we look at the unique solution 



(10) 



X{Ua) + A{x,v)Ua 
UA\d+(SM) = Id. 



0, (x,v) e SM 



The scattering data of the connection A is now the map Ca '■ d^{SM) U{n) defined 



as C A 



U 



A\d-{SM)- 



We can now rephrase Question I as follows: 



Question I. Let A and B be two unitary connections with Ca = Cb- Does there exist 
a smooth map U : M U{n) with U\9m = Id and B = U'^dU + V'^AUl 

Suppose Ca = Cb and define U := Ua{Ub)^^ '■ SM — )• U{n). One easily checks that 
U satisfies: 

' XU + AU -UB = 
U\d{SM) = Id. 

If we show that U is in fact smooth and it only depends on the base point x G M we 
would have an answer to Question I , since the equation above reduces to dU + AU — 
UB = and U\sm = Id which is exactly gauge equivalence. Showing that U only 
depends on x is not an easy task and it often is the crux of the matter in these type of 
problems. To tackle this issue we will rephrase the problem in terms of an attenuated 
ray transform. 

Consider W := U — Id : SM — )• C"^", where as before C"^" stands for the set of all 
n X n complex matrices. Clearly W satisfies 



(11) 
(12) 



XW + AW-WB 

W\d(SM) = 0. 



5 - A, 



We introduce a new connection A on the trivial bundle M x C"^" as follows: given a 
matrix R G C"^" we define A{R) := AR — RB. One easily checks that A is Hermitian 
if A and B are. Then equations (ITTl) and ([T^ are of the form: 

Xu + Au = -/, 

u\d{SM) = 0. 

where A is a unitary connection, / : SM — )■ is a smooth function linear in the 
velocities, u : SM — is a function that we would like to prove smooth and only 
dependent on a; G M and N = n x n. As we will see shortly this amounts to under- 
standing which functions / linear in the velocities are in the kernel of the attenuated 
ray transform of the connection A. 

First recall that in the scalar case, the attenuated ray transform laf of a function 
/ G C°°{SM,C) with attenuation coefficient a G C°°{SM,C) can be defined as the 
integral 



Iaf{x,v) 



t{x,v) 



f{(t>t{x,v))exp 



a{(f)s{x, v)) ds 



dt. 



{x,v) G d+{SM). 
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Alternatively, we may set laf := u\d+{SM) where u is the unique solution of the transport 
equation 

Xu + au = —f in SM, u\q_(^sm) = 0. 
The last definition generalizes without difficulty to the case of connections. Assume 
that A is a unitary connection and let / G C°°(S'M, C") be a vector valued function. 
Consider the following transport equation for u : SM C", 

Xu + Au = — / in SM, u\q_(^sm) = 0. 

On a fixed geodesic the transport equation becomes a linear ODE with zero initial 
condition, and therefore this equation has a unique solution denoted by u^ . 

Definition. The attenuated ray transform of f & C°°{SM,C^) is given by 

lAf ■= U^\d+{SM)- 

We note that Ia acting on sums of 0-forms and 1-forms always has a nontrivial 
kernel, since 

lA{dp + Ap) = for any p E C^{M, C") with p\9M = 0. 

Thus from the ray transform I^f one only expects to recover / up to an element having 
this form. 

The transform Ia also has an integral representation. Consider the unique matrix 
solution Ua '■ SM — )■ U {n) from above. Then it is easy to check that 

pt(x,v) 

lAf{x,v)= / U^\(j)t{x,v))f{(f)t{x,v))dt. 
Jo 

We are now in a position to state the next main question: 

Question II. (Kernel of Ia) Let (M, g) be a compact simple Riemannian manifold and 
let A be a unitary connection. Assume that / : SM — >■ C" is a smooth function of the 
form F{x) + aj{x)v^ , where F : M — )■ C" is a smooth function and a is a C"-valued 
1-form. If /a(/) = 0, is it true that F = and a = dAP = dp + Ap, where p : M — >■ C"" 
is a smooth function with p\dM = 0? 

As explained above a positive answer to Question II gives a positive answer to 
Question I. The next recent result provides a full answer to Question II in the two- 
dimensional case: 

Theorem 9.1. [33J Let M be a compact simple surface. Assume that f : SM — )■ is 
a smooth function of the form F{x) + aj{x)v^ , where F : M C"' is a smooth function 
and a is a C^-valued 1-form. Let also A : TM u{n) be a unitary connection. If 
^A^f) = 0, then F = and a = dAP, where p : M ^ is a smooth function with 

P\dM = 0. 

Let us explicitly state the positive answer to Question I in the case of surfaces: 

Theorem 9.2. [33j Assume M is a compact simple surface and let A and B be two 
unitary connections. Then Ca = Cb implies that there exists a smooth U : M ^U{n) 
such that U\dM = Id and B = U'^dU + U'^AU. 
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The proof of Theorem 19.11 is based on the ideas explained in Section [31 One intro- 
duces a suitable additional attenuation (twists with a positive line bundle) which adds 
positivity to the Pestov identity with a connection and then gauges the twist away via 
the key Theorem 13.11 

In the case of Euclidean space with the Euclidean metric the attenuated ray trans- 
form is the basis of the medical imaging technology of SPECT and has been extensively 
studied, see p!5] for a review. We remark that in connection with injectivity results 
for ray transforms, there is great interest in reconstruction procedures and inversion 
formulas. For the attenuated ray transform in with Euclidean metric and scalar 
attenuation function, an explicit inversion formula was proved by R. Novikov [2S]. A 
related formula also including 1-form attenuations appears in inversion formulas for 
matrix attenuations in Euclidean space are given in [HI and the case of hyperbolic 
space is considered in [S]. 

Various versions of Theorem 19.21 have been proved in the literature. Sharafutdinov 
[H] proves the theorem assuming that the connections are close to another con- 
nection with small curvature (but in any dimension). In the case of domains in the 
Euclidean plane the theorem was proved by Finch and Uhlmann [16] assuming that 
the connections have small curvature and by G. Eskin [H] in general. R. Novikov [29] 
considers the case of connections which are not compactly supported (but with suitable 
decay conditions at infinity) and establishes local uniqueness of the trivial connection 
and gives examples in which global uniqueness fails (existence of "ghosts" ) . 

For more on inverse problems for connections we refer to |30j . 

10. Anosov manifolds 

There are versions of the ideas in the previous sections in the context of closed 
manifolds. The first requirement is to have a notion that replaces the concept of 
simple manifold. It is easy to motivate this as follows. Simple manifolds have two 
characteristic properties: they have no conjugate points and and they are open in the 
C^-topology of metrics. Recall that a Riemannian manifold is said to have no conjugate 
points if any two points in the universal covering are joined by a unique geodesic. Hence 
it seems natural to seek an analogue by requiring that the metric is a C^-interior point 
among the set of all metric without conjugate points. 

Definition. A closed Riemannian manifold {M,g) is said to be Anosov if g belongs to 
the C'^-interior of the set of metrics without conjugate points. 

It turns out that the name "Anosov" is completely justified: (M, g) is Anosov if and 
only if the geodesic flow of g is Anosov in the sense of Dynamical Systems [IQ]. We 
will not give here the definition of an Anosov flow since it will not be explicitly needed 
and instead we refer the reader to [22] • 

From our definition it is clear that negatively curved manifolds are Anosov and that 
there are no Anosov metrics on tori since the only metrics without conjugate points on 
tori must be fiat [7]. 

The notion of "Jm is s-injective" makes sense for closed manifolds as follows: 
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Definition. We say that 1^ is s-injective if given any symmetric m-tensor f such that 

[ fMt),iit))dt = o 

Jo 

for every unit speed closed geodesic 7 : [0, T] — )■ M, then f is potential, i.e., there exists 
an (m — l)-symmetric tensor h such that f = dh. 

The tensor tomography problem for an Anosov manifold consists in proving that Im 
is s-injective for any m. There are numerous motivations for this, but perhaps the 
most notorious one is that of spectral rigidity which involves l2- In [19] Guillemin and 
Kazhdan proved that if (M, g) is an Anosov manifold such that I2 is s-injective then 
(M, g) is spectrally rigid. This means that if {gs) is a smooth family of Riemannian 
metrics on M for s G {—e,e) such that go = g and the spectra of — A^^ coincide up to 
multiplicity, 

Spec(-AgJ = Spec(-A3j, s G 
then there exists a family of diffeomorphisms tpg : M ^ M with ipo = Id and 

9s = Vs9o- 

Let us summarize what is known about the tensor tomography problem on an Anosov 
manifold. 

• /o and Ii are s-injective pT] : 

• I2 is s-injective for surfaces [M] ; 

• Im is s-injective for all m for non-positively curved manifolds [8j. 

11. Open problems 

In this section we mention some open problems related to tensor tomography. 

1. In the two dimensional case there is by now, as surveyed in this paper, a rather 
good understanding of the injectivity and range of the geodesic ray transform on 
tensor fields for simple manifolds. Important questions remaining are inversion 
formulas or reconstruction procedures of the solenoidal part of the tensor field 
from its geodesic ray transform. Certain inversion formulas were given in [38j . 
[23] for the case of constant curvature and close to constant curvature. 

2. In the case where dim (M ) > 3 it is not known whether is s-injective for a 
general simple manifold. This is known for Jq and Ji, but even the case of I2 is 
unknown at present. 

3. Support type theorems for the geodesic ray transform, where a tensor field is 
determined locally from its line integrals in a certain neighborhood, are known 
for the case of real analytic simple manifolds for J^, m = 0, 1, 2 [231 125]. Is it 
possible to extend these results to all simple manifolds? This has been done for 
Jo in three dimensions or higher |56j . 

4. The study of s-injectivity of the geodesic ray transform for non-simple manifolds 
is an important problem for which not much is known. Certain results are given 
in [43j, [45j. A microlocal analysis of Jq when the exponential map has fold 
type singularities was done in [M]. Injectivity, stability and reconstruction were 
proven for Jq in the case of three dimensions or higher when the manifold can be 
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foliated by strictly convex hyper surf aces [56]. This allows for conjugate points. 
The s-injectivity of I2 was analyzed in [SB] for a class of non-simple manifolds. 
However, the question if /q is injective on a compact non-trapping manifold 
with strictly convex boundary is open. 

5. The attenuated ray transform for an unitary connection on simple surfaces 
and Anosov surfaces has been extensively studied in [311 [33], HI]. It would be 
interesting to extend the results to the case of a non-unitary connection. 

6. For closed Anosov surfaces it is known that Im is s-injective for m = 0, 1, 2 . Is 
it true for all m? Also, is I2 s-injective for Anosov manifolds of dimension > 3? 

7. Finally, it would be natural to extend all this theory to more general classes of 
curves. By this we mean replacing geodesies by other natural set of curves like 
magnetic geodesies or geodesies of affine connections with torsion (thermostats). 
Concerning magnetic geodesies, the tensor tomography problem in 2D is solved 
in [1] using the ideas presented here and the results in [10]. See also [1]. 
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